Constraining the optical emission from the double pulsar system 

J0737-3039 



F. R. Ferraro^ R. P. Mi gnani^'^, C. Pallanca^, E. Dalessandro^, B. Lanzoni^, A. Pellizzoni^, 

; A. Possenti^ M. Burgay^ F. Camilo^ N. D'Amico^'^ A. G. Lyne^ M. Kramer^'^ R. N. 

Manchester^ 

O 

^ I ^ Dipartimento di Astronomia, Universitd degli Studi di Bologna, via Ranzani 1, 1-40127 
D ■ Bologna, Italy 

^ Mullard Space Science Laboratory, University College London, Holmbury St. Mary, 

Dorking, Surrey, RH5 6NT, UK 

I ^ Kepler Institute of Astronomy, University of Zielona Gora, Lubuska 2, 65-265, Zielona 

^ " Gora, Poland 

^ ^ INAF-Osservatorio Astronomico di Cagliari, localitd Poggio dei Pini, Strada 54, 1-09012 

2 ! Gapoterra, Italy 

ci . Columbia Astrophysics Laboratory, Columbia University, New York, NY 10027, USA 



6 



Dipartimento di Fisica, Universitd di Cagliari, Cittadella Universitaria, 1-09042 



^ . Monserrato, Italy 

^ ■ Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, The University 

(N 

(N 
O 



of Manchester, Manchester M13 9PL, UK 

^ MPI fuer Radioastronomie, Auf dem Huegel 69, 53121 Bonn, Germany 

^ CSIRO Astronomy and Space Science, Australia Telescope National Facility, PO Box 76, 

Epping NSW 1710, Australia 



c^: February 7, 2012 



ABSTRACT 



We present the first optical observations of the unique system J0737— 3039 
(composed of two pulsars, hereafter PSR-A and PSR-B). Ultra-deep optical ob- 
servations, performed with the High Resolution Camera of the Advanced Camera 
for Surveys on board the Hubble Space Telescope could not detect any optical 
emission from the system down to mF435w = 27.0 and mpeoew = 28.3. The esti- 
mated optical flux limits are used to constrain the three-component (two thermal 
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and one non-thermal) model recently proposed to reproduce the XMM-Newton 
X-ray spectrum. They suggest the presence of a break at low energies in the non- 
thermal power law component of PSR-A and are compatible with the expected 
black-body emission from the PSR-B surface. The corresponding efficiency of the 
optical emission from PSR-A's magnctosphcrc would be comparable to that of 
other Myr-old pulsars, thus suggesting that this parameter may not dramatically 
evolve over a time-scale of a few Myr. 

Subject headings: Stars:Binaries:General, Stars:imaging, Stars:Pulsar:Individual: 
PSR J0737-3039A, Stars:Pulsar:Individual:PSR J0737-3039B, techniques: pho- 
tometric 



1. INTRODUCTION 

The binary system J0737— 3039 is the first and (up to now) the only double-pulsar sys- 
tem ever detected (Burgay et al. 2003; Lyne et al. 2004). It is composed of PSR J0737-3039A 
(hereafter PSR-A) and PSR J0737-3039B (PSR-B), with spin periods of 22 ms and 2.7 s, 
respectively, orbiting each other with a 2.5 hr period along a mildly eccentric (~ 0.09) orbit. 

Radio observations are providing the timing data for exploiting this formidable system 
to test General Relativity (Kramer et al. 2006, Kramer & Stairs 2008). Moreover, they led 
to the discovery of a wealth of unprecedented phenomena, hke a short (~ 30 s) echpse of the 
radio signal from PSR-A (modulated at full or half period of PSR-B) when it passes at the 
superior conjunction (McLaughlin et al. 2004) and a strong modulation of the flux density 
from PSR-B along the orbit (Lyne et al. 2004). The secular variations of these peculiar 
effects have been used to put constraints on the geometry of the system and on its evolution 
due to General Relativity (Burgay et al. 2005; Breton et al. 2008, Perera et al. 2010). On 
the other hand, they clearly demonstrate interaction between PSR-B's magnetosphere and 
the spin-down luminosity released by PSR-A in the form of a relativistic wind of particles 
and Poynting flux In turn, this opens the possibility of probing the usually inaccessible 
magnetosphere and the local environment of a neutron star (e.g. Zhang & Loeb 2004; 
Lyutikov & Thompson 2005; Breton et al., in preparation). 

Following the original detection with Chandra (McLaughlin et al. 2004), J0737— 3039 
has been repeatedly re-observed in the X-ray band (Pellizzoni et al. 2004; Campana et al. 
2004; Kargaltsev et al. 2006; Chattcrjcc et al. 2007; Posscnti et al. 2008; Peflizzoni et al. 2008, 
hereafter P08). As for the timing, these observations led to the detection of X-ray pulsations 
from both pulsars: PSR-A shows a double-peaked proflle with a ~ 70% pulsed fraction along 
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the entire orbit (Chatterjee et al. 2007), whereas PSR-B displays pulsations with a complex 
profile for (at least) a quarter of the orbit around the ascending node (P08). As for the 
spectroscopy, the phased-averaged X-ray spectrum can be satisfactorily fit by two-component 
models (a power-law plus a black body model, or a double black body model; see P08 and 
Possenti et al. 2008, respectively), while the phase-resolved spectrum is most hkely described 
in terms of magnetospheric emission from PSR-A and thermal emission from both the polar 
caps of PSR-A and the bulk of the PSR-B surface, irradiated by its companion (P08). If 
confirmed, the latter model would represent an additional manifestation of the interaction 
between the PSR-A's spin-down flux and the magnetosphere of PSR-B. Moreover, the study 
of the cold black-body component would allow one to shed light on long-standing questions 
in the physics and electro-dynamics of neutron stars, like the thermal inertia time of the 
surface and the degree of anisotropy (and the geometry) of the emission of the spin-down 
luminosity. 

Given the unprecedented nature of the J0737— 3039 system and the peculiarities shown 
in both the radio and X-ray bands, deep observations of the pulsar system in the optical 
domain are needed to complete the study of its phenomenology. In particular, optical ob- 
servations of pulsars are crucial for constraining the tails of both thermal and non-thermal 
spectra observed in the soft X-ray regime and to unveil possible deviations (see Mignani 2011 
for a recent review). Indeed, spectral breaks or turn-offs are found to be quite common in 
the thermal/non-thermal optical-to-X-ray spectra of pulsars (e.g. Mignani et al. 2010) and 
are crucial both to study anisotropies in the temperature distribution on the surface of the 
neutron star and to track the different energy and density distribution of particles in distinct 
regions of the pulsar magnetosphere, where the non-thermal optical and X-ray radiation is 
produced. 

Here, we present the results of the first, deep, optical observations of the J0737— 3039 
system. The data reduction procedures are presented in §2J Results and discussion are given 
in ^ Summary and conclusions follow in ^ 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. Description of the observations 

The photometric data-set used for this work consists of deep, high-resolution images 
obtained with the High Resolution Channel (HRC) of the Advanced Camera for Surveys 
(ACS) on board the i/ubble Space Telescope (HST). The HRC consists of one detector 
with a pixel-scale of ~ 0'.'025 and a field of view (FOV) of ~ 29" x 26". The images have 
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been obtained on 2005 July 20-21 (Prop. ID 10411, P.I. Ferraro) in ACCUM mode with 
a standard multi-point dithering pattern and in two different bands: six images through 
the F606W fiher (A = 5887.53A; AA = 664.7A; roughly corresponding to the Johnson's V 
band) with an exposure time tgxp = 1250 s each, and eight images through the F435W filter 
(A = 4310. 97A; AA = 309. 5A; corresponding to the Johnson's B band) with an exposure 
time texp = 1300 s for six of them, and texp = 1250 sec for the remaining two. Three short 
exposures of 10 s each were also acquired in the F606W band. We retrieved the data from the 
HST science data archive, after on-the-fly re-processing by the latest version (v. 6. 1.0) of the 
ACS data reduction pipeline (calacs), which applies bias and dark subtraction, overscan 
and flat-field correction, and photometry calibration using the most updated calibration files. 

As shown in Figures [Hand [21 the pulsar system (with radio coordinates a = 07'^37™51!248 
and 6 = -30°39'40'.'714; Lyne et al. 2004) is located at only ~ 4" from a bright star 
{Bj = 13.9 according to the GSC2 catalogue; Lasker et al. 2008) also detected in the 2MASS 
survey (see Table [1]). In order to prevent it to severely contaminate the pulsar region, we 
used the HRC coronograph to mask this object during the long exposure acquisition. 

2.2. Astrometry 

Partially due to the unaccounted proper motions of the reference stars, the absolute 
accuracy of the HST astrometry can be worse than the nominal value (~ 0'.'3). Hence, in 
order to precisely locate the pulsar system in the ACS/HRC frame we applied the following 
approach. 

Since the brightest star is obscured by the coronograph in the long exposures, for the 
astrometric re-calibration of the HRC frame we used the three F606W short images. In par- 
ticular, we used the drizzled (drz) images produced by the ACS data reduction pipeline, which 
are already corrected for the well known geometric distorsion affecting the HRC (Pavlovsky 
et al. 2005). By using the IRAF tasks geomap and gregister we aligned the single images 
taking into account shifts and rotations and then, by using the iRAF task imcombine, we 
combined them and applied a cosmic-ray rejection algorithm. 

After a quick reduction aimed at just obtaining the instrumental coordinates of the stars, 
we cross-correlated our stellar list with the 2MASS point source catalogue (Skrutskie et al. 
2006). Only three 2MASS stars have been identified within the FOV of our HRC observations 
(see Table [1] and Figs. [T]and[2]), their positions (epoch 2000.0) being systematically offset 
with respect to the default HST astrometry by an average value of S{RA) = — 0'.'21 and 
6{Dec) = 0'.'17, of the same order of magnitude as found in previous works (e.g., Ferraro et 
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al. 2001). The image was thus registered on the 2MASS astrometric grid, which is tied to the 
international celestial reference frame to better than 0'.'15 (Skrutskie et al. 2006). Note that 
the brightest star (2MASS-1) in the FOV has a measured proper motion: fiacosS = — 42.4±4.8 
mas yr~^ and fig = +30.4±4.8 mas yr~^ in the UCAC2 catalogue (Zacharias et al. 2004), and 
fJ'acosS = — 32.5±4.1 mas yr~^ and fis = +29.9±4.1 mas yr~^ in the recently released PPMXL 
catalogue (Roeser et al. 2010). This proper motion would only yield a small displacement 
(of the order of O'.'IS - 0'.'2) between the epoch of the 2MASS positions (2000.0) and that 
of our observations (2005.55). Taking into account the rms of the coordinate difference 
and the uncertainty introduced by the proper motion of star 2MASS-1, we adopt 0'.'3 as a 
conservative estimate of the global position error of our astrometry. The resulting location 
of the pulsar system in the ACS/HRC frame is shown in Figs. [1] and |2) it turns out to 
be ^ 3'.'27 3'.'05 west and ^ I'.'IT south) away from the bright star 2MASS-1. We note 
that the pulsar proper motion, fiacosS = —3.82 ± 0.62 mas yr~^ fis = 2.13 ± 0.23 mas yr~^, 
only yields a negligible displacement during the ~ 2 year difference between the epoch of 
the reference radio position and that of our HST observations. 

2.3. Photometry 

In order to search for even the faintest sources, we performed the photometric analysis 
on the median images, by applying to the long exposures in each filter the same procedure 
described above. At the end of the procedure, no object could be detected by eye at the 
pulsar position, neither in the F606W nor in the F435W images. 

For a more reliable and deeper search of any possible counterpart, we estimated the ob- 
servational detection limit, i.e. the minimum counts over the background level that the PSF 
peak of the faintest star should have to be successfully detected by the threshold criterium. 
To this end, we first determined the local background level and its standard deviation {<7bkg) 
in a 1" X 1" region centered on the nominal pulsar system position. Indeed this is a critical 
task in this particular case, since the extended wings of star 2MASS-1 may affect the local 
background level, despite the use of the coronograph. Following a standard approach, the 
number of counts corresponding to 3 x abkg was adopted as the observational detection limit. 
We then performed a very accurate photometric analysis around the nominal position of the 
pulsar by using romafot (Buonanno et al. 1983), a semi-automatic package that models 
even complex PSFs through an analytic function (Moffat, Gaussian, etc.) combined with 
a numerical residual matrix, and that allows a visual inspection of the quality of the fit. 
In each filter, we first determined the shape of the PSF by modeling five bright and non 
saturated stars in the HRC FOV, by using a Moffat function combined with the numerical 
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residual matrix. We then used such a PSF to search for sources at more than 3 x ai,kg in 
a region centered on the pulsar system position. No signal was found above this limit even 
when applying a Gaussian smoothing to the images over a 3 x 3 pixel cell. We therefore 
conclude that J0737— 3039 is undetected in our ACS/HRC images. 

The magnitude of the faintest detectable object has been computed as the magnitude 
of a star having the adopted PSF and a number of counts at peak equal to 3 x (Jbkg- By 
adopting the Charge Transfer Efficiency correction (Maybhate et al. 2010) and the absolute 
flux calibration (from the updated PHOTFLAM coefficients available in the most recent 
SYNPHOT tables maintained at StScI) our 3abkg upper limits turned out to correspond to 
"^F435W = 27.0 and mpeoew = 28.3 in the STmag photometric system^ 

We finally transformed these magnitude limits into upper limits to the unabsorbed spec- 
tral flux. To do this, we flrst evaluated the extinction in the direction of the pulsar system, 
by using the relation of Predehl & Schmitt (1995) between the hydrogen column density iVn 
and the extinction coefficient Ay- While this relation is affected by uncertainties for close 
objects, due to the problems of modeling the interstellar medium at small distance from the 
Sun where microstructures weight more, we checked that, when using, e.g., the relation of 
Paresce (1984), the resulting extinction corrections are consistent within 0.01 magnitudes. 
Assuming A^'h = 6.9tli x lO^o cm'^ (P08)@ and Ay = 3.12 x E{B - V), we obtained a 
colour excess E{B — V) = 0.123. This was used to compute the extinction coefficients ap- 
propriate for the central wavelength of the adopted filters (Fitzpatrick 1999), thus obtaining 
^F435W = 0.52 and Apeoew = 0.35. Using the standard relation for the STmag system 
(magA = —2.5 log Fa — 21.1), we obtained the monochromatic fluxes -Ff435w = 9.3 x 10"^° 
erg cm^^ s^^ A~^ and -FVeoew = 2.4 x 10~^° erg cm~^ s^^ A^^. The corresponding quantities 
in units of the frequency have been obtained by considering the effective wavelength of the 
two pass-bands: Af435w = 4310.9702 A and Apeoew = 5887.536 A. Finally, according to 
standard unit conversions (e.g. Zombeck 2007), we derived the values of the unabsorbed 
spectral fluxes -Ff435w = 8.7 x 10~^ keV cm~^ keV~\ and -Fpeoew = 4.2 x 10~^ keV cm"^ 
s"^ keV-^ 



^As a sanity check, we used the empirical PSF to simulate two artificial stars of different magnitude 
(mpeoew — 28.5 and 28.3) directly in the images, at the nominal position of the pulsar system. By re- 
performing the entire photometric analysis on these images, we found that only the brighter star could be 
successfully detected and measured. 

^The quoted uncertainties in translates in a 0.08 mag uncertainty in Ay. 
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3. DISCUSSION 

The derived upper limits to the optical flux in the F435W and F606W fllters can be 
used to characterize the optical emission of the J0737— 3039 system (§3.3), as well as to 
constrain the spectral model used to flt the X-ray data (§3.2), which is also briefly reported 
below (§3.1) for the sake of clarity (a detailed discussion of the X-ray spectra of the system 
is given in P08). 

3.1. The X-ray spectrum 

By exploiting deep XMM-Newton observations performed on October 2006, P08 could 
perform phased-resolved spectroscopy of J0737— 3039 and found that a model with three 
components, a power law (PL) and two black-bodies (BBs), is required to flt both the phase- 
resolved spectra and the energy dependence of the pulsed fraction of PSR-A. According to 
this model, the PL (with photon index F = 3.3lo;2) is responsible for most of the pulsed 
emission originating in the magnetosphere of PSR-A (e.g. Zhang & Cheng 1999; Zhang 
& Harding 2000). The hotter and fainter BB (HBB, with kT = 134l{J eV) has a ~ 100 
m emitting radius [f| and it is likely associated with the polar caps of PSR-A, re-heated 
by back-flowing particles accelerated in its magnetosphere (e.g. Cheng & Ruderman 1980; 
Arons 1981; Zavlin & Pavlov 1998; Zavlin 2006; Bogdanov 2006). The colder and brighter 
BB (CBB, with kT = 32^4 eV) is associated with an emitting radius of ~ 15 km and 
plausibly arises from the bulk of the surface of one of the two neutron stars. The most likely 
explanation is that it is due to thermal radiation coming from the bulk of PSR-B's surface, 
after the conversion of ~ 2% of the spin-down luminosity of PSR-A intercepted by PSR-B 
at suitable orbital phases (see e.g. Zhang & Loeb 2004; Lyutikov & Thompson 2005). This 
is also consistent with the observed soft pulsed flux from PSR-B (P08). 



■^The quoted radii of the BB components have been obtained by P08 under the assumption of a distance 
of 0.5 kpc, inferred from the dispersion measure of the two pulsars (Lyne et aL 2004) and a model for the 
distribution of electrons in the Galaxy (Cordes & Lazio 2002). We used this distance as a reference value. 
Within the uncertainties in spectral modeling and fitting, the resulting parameters are compatible with those 
derived from the still preliminary determinations of the annual geometric parallax of the system (0.2-1 kpc, 
Kramer et al 2006; - 1 kpc, Deller et al. 2009). 
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3.2. Optical vs. X-rays 

In order to use the derived optical flux upper limits for constraining the spectral model 
discussed above and getting new clues on the physics of the system, it is necessary to extrap- 
olate the PL and BB components into the optical domain. In agreement with the evidence 
found in the optical spectra of a few Myr old pulsars (Mignani 2011), a contribution in the 
optical band from the magnetospheric emission can be expected. Since PSR-A's spin-down 
energy is ~ 3000 times higher than PSR-B's one {Ea = 5.9 x 10^^ erg s'^; Eb = 1-7 x 10^° 
erg s^^; Lyne et al. 2004), we assume that the former would dominate any magnetospheric 
optical emission from the system, as it does in the X-rays. On the other hand, we note that 
the actual broad-band spectrum of PSR-A is, in principle, compatible with comptonised BB 
models (see e.g. Nishimura et al. 1986), which are not expected to provide a non-thermal con- 
tribution at optical wavelengths. While further X-ray observations are required to improve 
the photon count statistics and constrain complex multi-component phase-resolved spectral 
models (including comptonisation scenarios), here we conservatively take into account the 
possibility of a contribution of the non-thermal component even at low energies. 

The extrapolation of the PL and BB components into the optical domain is shown in 
Fig. [3l together with the optical flux upper limits computed in §2.31 While the expected 
contribution of the HBB is at least four orders of magnitude lower than the derived opti- 
cal values, the extrapolation of the CBB component only slightly undershoots the optical 
flux upper limits, the difference being smaller than a factor of 10 after accounting for the 
90% confldence level uncertainty of the spectral parameters. Finally, the extrapolated non- 
thermal PL component of the X-ray spectrum overshoots the optical flux upper limits by 
about six orders of magnitude. 

As a flrst consideration, our optical upper limits indicate that either the PSR-A mag- 
netosphere does not provide a non-thermal contribution at low energies (as expected in the 
comptonised BB model), or a break at these wavelengths is present jj Indeed, such breaks 
are quite common in the optical-to-X-ray magnetospheric spectra of pulsars above a wide 
age range (e.g., Mignani et al. 2010) and are probably indicators of a complex particle den- 
sity and velocity distribution in the neutron star magnetosphere. However, in order to be 
compatible with the expected optical emission from the CBB component (see Fig. [3]), the 



*Note that if the X-ray PL component was extrapolated to optical wavelengths with the same index, it 
would imply an emitted power in the optical band which is a significant fraction (0.1-0.2) of the measured 
PSR-A's spin-down luminosity, in striking contrast with what observed so far in the pulsar population (see 
§3.3). This argument further supports the presence of a spectral break between the optical and the X-ray 
domains. 
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expected spectral break in PSR-A's spectrum has to be very sharp if it occurs at energies 
lower than ~ 0.1 keV. 

In the case of a null or negligible optical contribution of the PL component, the fact 
that the measured upper limits are consistent with the expectations from the combined BB 
components represents a partial validation of the model and the adopted spectral parameters. 

By neglecting the contribution of PSR-A's polar caps (i.e., the HBB component), the 
comparison between the optical upper limits and the expected CBB contribution more quan- 
titatively suggests that the PSR-A's magnetospheric optical emission cannot exceed the sur- 
face thermal emission from PSR-B by a factor larger than ~ 5. Conversely, by assuming 
that the PL spectral break is such to make any magnetospheric optical emission negligible 
with respect to the thermal one, we can derive a qualitative upper limit to the temperature 
of the bulk of PSR-A surface (i.e. away from the hot polar caps) to account for the possible 
flux excess with respect to the expected CBB contribution. For instance, in order to be 
compatible with the non-detection of a third BB component (cooler than the CBB) in the 
observed X-ray spectrum, for an emitting radius of 15 km, the bulk of PSR-A surface should 
be at a temperature kT < 20 eV. This is consistent with its spin-down age (~ 200 Myr) and 
a passive cooling scenario. 

Spectral data in the ultraviolet would be useful to constrain the presence and energy 
of a break at low energy in the PL spectrum of PSR-A, as well as the temperature of the 
CBB component produced by PSR-B. Unfortunately, J0737— 3039 has not been observed by 
GALEX (Martin et al. 2005) and the XMM-A^ew;ton Optical/UV monitor camera (Mason 
et al. 2001) has been used only in the V band, while the relatively recent discovery of the 
system prevented its observation within the pulsar survey carried out by EUVE (Korpela & 
Bowyer 1998). 

3.3. Optical emission efficiency 

The F606W upper limit corresponds to an extinction-corrected optical luminosity for 
the pulsar system //F606W ^ 4.8 x lO^^figQQ erg s~^ (where d^oo is its distance in units of 500 
pc), with a < 7% uncertainty due to the error on the extinction correction computed from 
A'h (see §2.31) . In the limit case that PSR-A's emission is mostly magnetospheric, and under 
the assumption that its contribution to I^Feoew is at most a factor of 5 larger than that 
of the PSR-B's thermal emission (see above), we can derive the optical emission efficiency: 
Vopt,A = Lfqogw/Ea < 8 X 10"*^ (igQQ. 

Such a value of rjopt implies that the optical emission efficiency of PSR-A might still be 
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compatible with that of most pulsars with strong magnetospheric emission (e.g., Zharikov 
et al. 2006), accounting for the uncertainties on the distance and for the difference in the 
reference spectral banco. This might indicate that the optical emission efficiency does not 
dramatically evolve with the pulsar's spin-down age and remain relatively high for pulsars as 
old as ~ 200 Myr, like PSR-A (Lorimer et al.2007). Unfortunately, given the small number 
of pulsars identified in the optical domain, of which only a handful are older than 1 Myr 
(see Mignani 2011), it is very difficult to outline a clear evolutionary pattern. For instance, 
for the 3 Myr old PSR B1929+10 (Mignani et al. 2002) the optical luminosity, z/ entirely 
magnetospheric, would imply rjopt ~ 4.7 x 10^^, i.e. above the upper limit derived for PSR- 
A. Under the same assumption, for the ~ 4 Myr old PSR B1133+16 (Zharikov et al. 2008) 
and for the 17 Myr old PSR B0950+08 (Zharikov et al. 2004) we derive r]opt ~ 10"^ and 
~ 1.3 X 10"^ respectively, while for the ~ 166 Myr old PSR J0108-1431 (Mignani et al. 
2008; 2011), whose optical identification is however unconfirmed, rjopt would be as high as 
^ 5 X 10"^ 

Surprisingly enough, these values are comparable to those of much younger and brighter 
pulsars, like the Crab and PSR B0540— 69. However, we warn here that, while the optical 
emission of young pulsars is entirely magnetospheric (e.g., Mignani et al. 2010), that of the 
old pulsars partially (if not entirely) arises from the neutron star surface. Thus, the above 
quoted values of riopt are rather uncertain, and optical spectra would be needed to disentangle 
the contributions of the magnetospheric and thermal components. Indeed, in the case of PSR 
J0108— 1431, whose spin-down age is comparable to that of PSR-A, the optical spectrum of 
its (candidate) counterpart is compatible with the Rayleigh- Jeans tail of a kT ~ 0.43 d^^Q eV 
BB for a 13 km neutron star radius (Mignani et al. 2008; 2011), implying that the bulk of the 
optical emission arises from the neutron star surface and that the actual emission efficiency 
is much lower than 5 x 10~^. On the other hand, since PSR-A has a ~ 1000 times larger 
spin-down energy, its magnetospheric emission is probably much stronger than that of PSR 
J0108— 1431 and, possibly, it is dominant over the thermal component. Thus, the derived 
luminosity upper limits do not rule out that the (yet undetected) optical emission from the 
pulsar system could be of magnetospheric origin and almost entirely associated with PSR-A. 

4. Summary and conclusions 

We performed the first ever follow-up optical observations of the J0737— 3039 system. 
Our deep ACS/HRC observations could not detect optical emission from the system down 



'''The efficiencies quoted in Zharikov et ah (2006) are computed in the B band. 
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to limiting magnitudes of mF435w = 27.0 and mpeoew = 28.3, corresponding to upper limits 
to the unabsorbcd spectral fluxes -Ff435W = 8.7 x 10~^ keV cm~^ s~^ keV~^ and i^reoew = 
4.2 X 10^^ keV cm"^ s"^ keV^^. In the framework of the PL+2BB model that best fits the 
XMM- Newton spectrum (P08), the derived flux upper limits suggest a break at low energies 
in the PL component, similar to what observed in other pulsars (e.g., Mignani et al. 2010), 
and are consistent with the BB temperature and emitting area of PSR-B. The upper limit 
to the optical luminosity of the pulsar system (Lpeoew ^ 4.8 x 10^^d|oo ^^S ^~^) imphes 
that the optical emission efficiency of PSR-A, which has the highest spin-down power, would 
be comparable to that of Myr-old pulsars, suggesting that their optical magnetospheric 
emission does not dramatically evolve over a time-scale of a few Myr. Deep observations in 
the ultraviolet with the HST and its refurbished instrument suite will probably offer higher 
chances of detection and, thus, of better constraining the thermal spectrum of PSR-B and 
the magnetospheric spectrum of PSR-A. 
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Table 1: Absolute positions of the three 2MASS stars (Skrutskie et al. 2006) used to compute 
the astrometric solution in the ACS/HRC frame (§2.2). 



NAME 


R.A. 


Dec. 




(J2000) 


(J2000) 


2MASS-1 


07 37 51.48 


-30 39 39.57 


2MASS-2 


07 37 51.85 


-30 39 43.41 


2MASS-3 


07 37 52.51 


-30 39 54.05 




2MASS-1 

® J0737-3039 



2MASS-2 _ 



2MASS-3 I 1 

m 

Fig. 1. — ACS/HRC image of the region around J0737— 3039 taken through the F606W filter 
(7500 s) after astrometric re-calibrations (see §2.21) . The small, white squares superimposed 
to the three brightest stars in the field correspond to the 2MASS coordinates given in Table 
[U The cross corresponds to the estimated location of the J0737— 3039 system, and the circle 
marks the fiducial uncertainty (0'.'3) on our astrometric solution. 




Fig. 2. — Zoomed region around the location of J0737— 3039 in the ACS/HRC image ref- 
erence frame. The structures observed at the position of the bright star 2MASS-1 are the 
effects of the HRC coronograph. 
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Photon Energy [keV] 



Fig. 3. — Extrapolation into the optical regime of the model that best fits the XMM-Newton 
data of the pulsar system: a cooler BB and a hotter BE plus a PL. The solid lines are the 
best fit model components, while the dot-dashed lines correspond to the 90% confidence level 
uncertainties of the best fit. The grey region marks the location and uncertainties of the 
X-ray measurements. The two arrows indicate the F435W and F606W-band upper limits 
obtained from the HST data. Since the uncertainty on the extinction correction derived 
from the assumed A^h implies only a difference of < 7% in flux, we have neglected it when 
plotting the upper limit values. 



